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We use the nonperturbative linear 8 expansion method to evaluate analytically the coefficients 
ci and c'2 which appear in the expansion for the transition temperature for a dilute, homogeneous, 
three dimensional Bose gas given by T c = T {l + cicm 1/3 + [c 2 ln(an 1/3 ) + C2]a 2 n 2/3 + 0(o 3 n)}, where 
To is the result for an ideal gas, a is the s-wave scattering length and n is the number density. In a 
previous work the same method has been used to evaluate ci to order-5 2 with the result ci = 3.06. 
Here, we push the calculation to the next two orders obtaining c 1 — 2.45 at order-J 3 and C\ = 1.48 
at order-5 4 . Analysing the topology of the graphs involved we discuss how our results relate to 
other nonperturbative analytical methods such as the self-consistent resummation and the 1/N 
approximations. At the same orders we obtain d{ = 101.4, c'J, — 98.2 and c% = 82.9. Our analytical 
results seem to support the recent Monte Carlo estimates ci = 1.32 ± 0.02 and d£ — 75.7 ± 0.4. 

PACS numbers: 03.75.Fi, 05.30.Jp, 11.10. Wx 



I. INTRODUCTION 

Recently, the evaluation of the critical temperature for interacting dilute homogeneous Bose gases has been the 
interest of many theoretical works. For this purpose, the starting model is the one used in the analysis of a gas of 
interacting boson particles, described by a complex scalar field tp, with a local interaction characterized by the s-wave 
scattering length a and Euclidean action which, in natural unities, can be written as 



J\t J d 3 x{v>*(x,r)(|--^V 2 )^x,T)-^ . (1.1) 



The field vb can be decomposed into imaginary-time frequency modes ipj (x, uij ) , with discrete Matsubara frequencies 
LUj = 2-KjJB and j being an integer whereas (3 is the inverse of the temperature. At the early stages of solving this 
problem jl| the non-zero Matsubara frequency modes have been integrated out generating a reduced three dimensional 
0(2) scalar theory. This procedure was justified on the grounds that near the transition the non-zero Matsubara modes 
decouple and one is left with an effective action given by 



S 3d = f3j d 3 xU*o (-^ v2 -^o + ^[^o] 2 | • (1-2) 

'X: 

Despite this simplification the problem remains non-trivial since ordinary perturbation theory cannot be used to treat 
the model at the phase transition due to the severe infrared divergences for the zero frequency modes ^>o at the critical 
point, originating the breakdown of conventional perturbation theory. Different nonperturbative methods, some of 
which are currently used in quantum field theories, have then been used to compute the transition temperature. The 
analytical methods include the self-consistent resummation (SCR) used by the authors of Ref. jlj, the 1/N expansion 
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used at leading order (l/N-hO) by Baym, Blaizot and Zinn-Justin and at next to leading order (1/iV-NLO) by 
Arnold and Tomasik Q as well as the linear (5-expansion (LDE) employed by some of the present authors in Ref. 
P). The numerical methods used mainly Monte Carlo lattice simulations (MCLS) like the ones employed recently by 
Arnold and Moore § and by Karshunikov, Prokof'ev and Svistunov ||. Most of those calculations predicted that, in 
the dilute limit, the shift of the critical temperature of the interacting gas, T c , as compared to the critical temperature 
for an ideal gas, To, AT C = T c — To, behaves as 

^^c.an^ + oUn 2 / 3 ) , (1.3) 

where n is the number density, c\ is a numerical constant and the critical temperature for an ideal gas is given as 
usual by 



2tt 

To = — 

TO 



C(3/2) 



(1.4) 



The constant C\ in Eq. (1.3) is directly related to the contributions from the zero mode Matsubara frequencies 
and therefore can only be computed from nonperturbative methods. Some recent numerical applications predicted 
values for c\ which are close to 1.30 (MCLS, ^ ||). On the other hand, the analytical applications mentioned above 
predicted the values 2.90 (SCR, §), 2.33 (1/AT-LO, §), 1.71 (1/JV-NLO, §]) and 3.06 (LDE, @). Additionally, the 
authors of Ref. |?J have also argued that a logarithmic term appears at order-a 2 in Eq. ( |l.3| ). They have shown that 
this term is of the form c' 2 a 2 n 2 / 3 ln(an 1//3 ) and also estimated, using large- N arguments, the value of the numerical 
coefficient c' 2 - Recently, Arnold, Moore and Tomasik Q have argued that when naively going from the original action 
(Se) to the reduced action (53d) by ignoring the effects of non-zero frequency modes one misses the effects that 
short-distances and/or high-frequency modes have on long-distance physics. For T c (n) at second order these effects 
can be absorbed into a modification of the strengths of the relevant interactions which means that one should consider 



the more general form for the reduced effective action Eq. (1.2) 



Srfffofo,^] =[3jd 3 X {Vo (~ Z 0^ V2 " W») 4>0 + lr M 2 + O [^o|V^| 2 , (V>*V>) 3 ] } + /^vacuum , (1.5) 

where Z^ is the wave-function normalization function, /13 incorporates the mass renormalization function, Z a incor- 
porates the vertex renormalization function and T vacuum represents the vacuum energy contributions coming from the 
integration over the nonstatic Matsubara modes. The O [^q -0ojV-0o| 2 : (V'o "^o) 3 ] terms represent higher order inter- 
actions in the zero modes of the fields. As emphasized in Ref. B, these terms will give contributions to the density 
of order a 3 and higher and therefore do not enter in the order-a 2 calculations. By matching perturbative order-a 2 
results obtained with the original action Se and the general effective action S e g, the authors of Ref. || were able to 
show that the transition temperature for a dilute, homogeneous, three dimensional Bose gas can be expressed at next 
to leading order as 



AT C 1/3 
= aan ' 



c^lmW/ 3 ) + 4' a 2 n 2 ' 3 + O (a 3 n) . (1.6) 



A similar structure is also discussed in Ref. g|. As far the numerical coefficients are concerned, the exact value for 
4, c ' 2 = -64ttC(1/2)C(3/2)- 5 / 3 /3 ~ 19.7518, was obtained using perturbation theory ||. The other two coefficients 
cannot be obtained perturbatively but they can, through the matching calculation, be expressed in terms of the two 
nonperturbative quantities n and 1Z which are, respectively, related to the number density (ipQifio) an d to the critical 
chemical potential /x c , as shown below. The actual relation in between the two nonperturbative coefficients and these 
physical quantities is given by S 



and 



d = -128tt 3 [C(3/2)]- 4 / 3 K , (1.7) 



|[C(3/2)r 5 / 3 t4' + I[C(3/2)]- 8 / 3 (192^ 3 K ) 2 + ^C(l/2)[C(3/2)]- 5 / 3 lnC(3/2) , (1.8) 
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where 6' 2 ' in Eq. (1.8) is given by 



&2 = 32tt 



ln(1287r 3 ) 



72tt 2 TZ - 96tt 2 k 



C(V2) 



T7Z [C(1/2)]2 



(1.9) 



with if 2 = —0.13508335373. The quantities k and 1Z are related to the zero Matsubara modes only. Therefore, they 
can be nonperturbatively computed directly from the reduced action S c s which, as discussed in the numerous previous 
applications, can be written as 



5*0 



d A x 



i|V0| 2 + ir barc 2 + J(0 2 ) 2 



(1.10) 



where (j) = (<f>i,<f>2) is related to the original real components of ipo by ?/>o( x ) = [</>i(x) + ifofa)], ?*bare — 

—2m\i-sl and u = A&itamT(Z a / Z^ 2 ). The vacuum contribution appearing in Eq. ( |l.5j ) will not enter in the specific 
calculation we do here. 



The three dimensional effective theory described by Eq. (1.1C) is super renormalizable 1 requiring only a mass 
counterterm to eliminate any ultraviolet divergence. In terms of Eq. ( 1. 10| ) , the quantities k and 1Z appearing in Eqs. 
( Q - (|T|) are defined by g 



(1.11) 



and 



K: 



S(0) 



(1.12) 



where the subscripts u and in Eq. ( [Flip mean that the density is to be evaluated in the presence of interactions 
and in the absence of interactions, respectively, and 2(0) is the self-energy with zero external momentum. Since they 
dependent on the zero modes their evaluation is valid, at the critical point, only when done in a nonperturbative 
fashion. As discussed in the next section the relation between r c and S(0) comes from the Hugenholtz-Pines theorem 
at the c ritic al point. 

Eq. (|L^) is a general order-a 2 result with coefficients that, therefore, depend on nonperturbative physics via k 
and 1Z. In principle, to evaluate these two quantities one may start from the effective three-dimensional theory, given 
by Eq. (1.10), and then employ any nonperturbative analytical or numerical technique. In general, the analytical 
nonperturbative methods give a prescription so as to select and sum an infinite number of contributions belonging 
to a given class. For example, the infinite subset that contains only direct (tadpole) contributions represents the 
Hartree approximation, whereas exchange contributions are also taken into account in the Hartree-Fock approximation. 
In practice the sum is achieved by using a modified ("dressed") propagator to evaluate physical quantities. The 
nonperturbative results are then generated by solving self-consistent equations. However, in resumming calculations 
the bookkeeping and renormalization may become a problem beyond leading orders. 

Another popular analytical nonperturbative technique is the 1/N expansion jl^, |ll[ where one sums infinite subsets 
of contributions whose order is labeled by 0(l/N n ) where N is the number of field components. In general, the leading 
order contribution is easily evaluated and may reveal interesting nonperturbative physics, at least from a qualitative 
point of view, apart for providing an "exact" result within the large-./V limit. A nice illustration is provided by its 
application, for example, to the Gross-Neveu model at zero temperature, where the issue of chiral symmetry breaking 
as well as asymptotic freedom were investigated Jl^] . From a quantitative point of view the leading order may not be 
sufficient and lead to errors since N is finite and not too large in most cases. An example of this case is illustrated by 
treating the same Gross-Neveu model at finite temperature, where the leading order large- N calculation predicts a 
finite value for the critical temperature at which chiral symmetry restoration takes place, in contradiction to Landau's 
theorem for phase transitions in one space dimension (l3|. 

In practice, going to higher orders can be a difficult task. Nevertheless, the 1/N ranks as a good method to investi- 
gate nonperturbative physics as shown in many applications. In particular, the results provided by this approximation 



1 Recall that the coupling constant, u, has dimensions of mass in natural unities. 



4 



for the interacting Bose gas case, where N = 2, are surprisingly good already at leading order M. Good numerical 
results can also be obtained with self-consistent methods despite some potential problems as discussed in Ref. [||. 
The numerical calculations use mainly Monte Carlo lattice techniques and many different results, for the interacting 
Bose gas critical temperature problem, were generated in this way. The differences arise mainly from the way the 
theory is put on the lattice, the size of the lattice, the way the continuum limit is taken and other issues. As already 
mentioned, two recent works seem to have settled this question [|[ ||. 

Here we shall present, and then apply, an alternative analytical nonperturbative method known as the linear 5 
expansion (LDE) fl4|, [l5[| (for earlier works, see for instance, Ref. JT6|), which is closely related to the variational 
perturbation theory Jl7[ and the Gaussian effective potential p8| . This same method re-appeared under the name 
of optimized perturbation theory Jl9| . The main attractive feature of this approximation is the fact that the actual 
evaluation of a physical quantity, including the selection of the finite subset of relevant contributions at each order, is 
done exactly as in perturbation theory. It is then easy to control and explicitly evaluate one by one each of the reduced 
number of contributions appearing at each order. The implementation of the renormalization procedure follows the 
one performed in most quantum field theory textbooks [^0| . After the usual perturbative manipulation one generates 
nonperturbative results through an optimization procedure, as we will discuss in the next section. 

This work is organized as follows. In Sec. II we present the method and illustrate it with a simple application 
to the pure anharmo nic os cillator. In the same section we implement the met hod in the e ffecti ve three dimensional 
theory given by Eq. (1.10) in order to evaluate the constants K and 1Z, Eqs. (1.11) and (1.12). The quantity r c is 
then evaluated in Sec. Ill whereas ((f) 2 ) is evaluated in Sec. IV. The optimization procedure is carried out in Sec. V 
where the numerical results are presented and compared with some of the recent results. We present our conclusions 
in Sec. VI. All contributions, which include difficult five-loop Feynman diagrams with arbitrary N, are explicitly 
evaluated by brute force without recurring to any approximations. An appendix is included to show the details of the 
calculations of these higher order terms. To our knowledge, some of them have not been evaluated in this way before. 



II. THE METHOD AND ITS APPLICATION TO THE INTERACTING BOSE GAS PROBLEM 



A. The linear 5 expansion 



The linear S expansion (LDE) was conceived to treat nonperturbative physics while staying within the familiar 
calculational framework provided by perturbation theory. In practice, this can be achieved as follows. Starting from 
an action S one performs the following interpolation 



S ^ S s = SS + (1 - 6)S (r]) , (2.1) 

which reminds the trick consisting of adding and subtracting a mass term to the original action. One can readily 
see that at 8 = 1 the original theory is retrieved. This parameter is really just a bookkeeping parameter and some 
authors do not even bother considering it explicitly as we do pl| . The important modification is encoded in the 
field dependent quadratic term Sq(t}) that, for dimensional reasons, must include terms with mass dimensions (77). In 
principle, one is free so as to choose these mass terms and within the Hartree approximation they are replaced by 
a direct (or tadpole) type of self-energy before one performs any calculation. In the LDE they are taken as being 
completely arbitrary mass parameters which will be fixed at the very end of a particular evaluation. One then formally 
pretends that 5 labels interactions so that So is absorbed in the propagator whereas 6S0 is regarded as a quadratic 
interaction. So, one sees that the physical essence of the method is the traditional dressing of the propagator to be 
used in the evaluation of physical quantities very much as in the Hartree case. What is different in between the two 
methods is that with the LDE the propagator is completely arbitrary while it is constrained to cope only with direct 
terms within the Hartree approximation. So, within the latter approximation the relevant contributions are selected 
according to their topology from the start. 

Within the LDE one calculates in powers of S as if it was small. In this aspect the LDE resembles the large- A 
calculation since both methods use a bookkeeping parameter which is not a physical parameter like the original 
coupling constants and within each method one performs the calculations formally working as if N — > 00 or 6 — > 0, 
respectively. Finally, in both cases the bookkeeping parameters are set to their original values at the end which, in 
our case, means 6 = 1. However, quantities evaluated with the LDE dressed propagator will depend on rj unless one 
could perform a calculation to all orders. Up to this stage the results remain strictly perturbative and very similar 
to the ones which would be obtained via a true perturbative calculation. It is now that the freedom in fixing 77 
generates nonperturbative results. Since rj does not belong to the original theory one requires that a physical quantity 
$ calculated with the LDE be evaluated at the point where it is less sensitive to this parameter. This criterion, known 
as the Principle of Minimal Sensitivity (PMS), translates into the variational relation |L8) 
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(2.2) 



The optimum value fj which satisfies Eq. ( |2.2|) must be a function of the original parameters including the couplings, 
which generates the nonperturbative results. The conver gen ce p roperties of this method has been rigorously proved 
in the context of the anharmonic oscillator (AO) |2l], p2[ p3|, p4|j . Very recently, Kneur and Reynaud |2j| claimed to 
have proved the convergence of this method in rcnormalizablc quantum field theories. These are very encouraging 
results for the present application which uses a renormalizablc effective model which shares many similarities with 
the pure AO. Let us quickly illustrate how this method works by considering the anharmonic oscillator described, in 
Minkowski space, by 



(2.3) 



If one sets m — in the relation above the model describes the pure anharmonic oscillator which cannot be treated 
by usual perturbation theory. Let us first consider the ground state energ y d ensity whose exact result, £ oxact = 
A 1 / 3 0. 420804974478 . . ., has been calculated by Bender, Olaussen and Wang p6[ . Following Eq. (2A) one may write 
the interpolated action as 



A = i(a o 0) 2 -^V-4^ 4 + 4 7? ^ 2 



from which one obtains the perturbative order-(5 result |2l| 



(2.4) 



+ OC 



Z7T Z 



+ °° dp y 2 5} 3 
2?t p 2 — r\ 2 4 



dp 



1 



2tt p 2 — rf 



0{8 2 ) 



(2.5) 



Now, setting 8 — 1 and applying the PMS optimization procedure one gets 



which is a self-consistent mass gap equation. It can be easily checked that with this solution one resumms exactly 
the same contributions that would appear in the usual Hartree approximation. The same procedure will capture the 
physics which arises from exchange terms at order-5 2 where the first contribution of this type appears together with 
order-5 2 direct (Hartree) contributions. Moreover, as shown in other applications [^tJ, the result furnished by Eq. 
( |2.6|) remains valid at second order if one considers only the direct terms and this pattern is valid at any order in 
8. The actual value predicted at this lowest order is £W — ^Hartree ^ A 1 / 3 0.429 which is only about 2% greater 
than the exact result. As shown in Ref. (2l[] this result can still be improved as one goes to higher orders. Here, we 
shall be mainly concerned with the nonperturbative evaluation of the vacuum expectation value (4> 2 ). This quantity, 
whose exact result is (0 2 ) cxact = A" 1 / 3 0.456119955748 . . . §§, was also evaluated in Ref. The optimum values 

were obtained with fj values coming from its direct optimization and also from the optimization of £ . At order-(5 the 
value (0 2 ) (1) = A~ 1 / 3 0.446456 was obtained from the direct optimization (fj = 1.259921) and (</> 2 ) (1) = A~ 1 ' /3 0.436789 
was obtained from the injection of fj — 1.14471, which was generated by the optimization of £. One then sees that 
the optimum (<j> 2 ) numerical values generated by the two optimization procedures are very similar which could be 
expected since, at each order, the diagrams which contribute to {<p 2 ) and £ have the same structure. 

At this stage it should be clear how nonperturbative results may be generated, through the variational PMS 
procedure, from the perturbative evaluation of physical quantities. As already mentioned, the effective model to be 
considered in the sequel for the description of the dilute Bose gas temperature bears may similarities with the AO. The 
main differences being the number of space-time dimensions concerning each case (which means that one has to deal 
with ultraviolet divergences in three dimensions) and the fact that the former is used to investigate a phase transition. 
Technically, as we shall see, this translates into extra difficulties due to the Hugcnholtz-Pincs theorem which washes 
out direct (tadpole) contributions, meaning that the first non-trivial contributions to (cf) 2 ) start at the three-loop level 
via two-loop self-energies. Apart from the quantum mechanical applications [|l[ |22], |23| , p4j , the LDE was successfully 
applied to the description of mesoscopic systems p9[, nuclear matter properties J27], phase transitions in the scalar 
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A</> 4 model |3C], [u| as well as in the Gross-Neveu model p3| , investigation of chiral symmetry phenomena in QCD |3j| 
and in the determination of the equation of state for the Ising model |34| . It is worth mentioning that the application 
of the LDE to the scalar 0(N) x O(N) model |35) has allowed to investigate the nonperturbative phenomenon of 
symmetry nonrestoration at high temperatures further than it was possible with other standard nonperturbative 
methods. 

The first application of this method to the present problem was performed in Ref. where only the first non- 
trivial contribution, which appears at order-<5 2 , was considered. A successful extension to the ultra-relativistic case 
was performed by Bedingham and Evans in Ref. |36| . 



B. The interpolated theory for the zero frequency Matsubara modes 



One can now write the interpolated version of the effective model described by Eq. (1.10). Before doing that let 
us rewrite rbarc = r + A where A is a mass counterterm coefficient. This countertcrm is the only one effectively 
needed within the modified minimal subtraction (MS) renormalization scheme which we will adopt here. Then, one 
can choose 



So = \[m 2 +V 2 cl> 2 ] , (2.7) 



obtaining 



S s = 



d 3 x 



\m 2 + \v 2 <? 



Su . 2 2 8 

—r(4>) + ~ 
4! v ' 2 



(2.8) 



Note that we have treated r (r c at the critical point) as an interaction, since this quantity has a critical value which is 
at least of order 8. The Feynman rules for this theory, in Euclidean space, are — 8r, 8rj 2 and — 8 As, for the quadratic 
vertices and — Su for the quartic vertex. The propagator is given by 



G (0) (p) = [p 2 +»7 2 r 1 . (2.9) 

The corresponding diagrams for these rules are shown in Fig. 1. Note that rj acts naturally as an infrared cutoff so 
we do not have to worry about these type of divergences. By introducing only quadratic terms the LDE interpolation 
does not alter the polynomial structure, and hence the renormalizability, of the theory. 

In general, the counterterm coefficients appearing in the interpolated theory have a trivial dependence on the 
bookkeeping parameter and the renormalization process can be consistently achieved with the interpolated theory 
exactly as in ordinary perturbation theory. Once inserted into a diagram, the extra quadratic vertex proportional to 
8rf brings in more propagators decreasing the ultraviolet degree of divergence. We point out that renormalization 
should be carried out before the optimization process to ensure that the optimum value fj is a finite quantity. The 
interested reader is referred to Refs. ||(], |35| for more details concerning renormalization within the LDE. 

Requiring that at the critical temperature the original system must exhibit infinite correlation length, means that, 
at T c and 8 = 1 (the original theory) , the full propagator (p) , given by 



G^(p)=[p 2 + r 1 2 + 8r-8r 1 2 + Eill(p) 
must satisfy G^(0) -1 = 0, which implies 



(2.10) 



JrW = -Egl (0) . (2.11) 



The above equation is equivalent to the Hugenholtz-Pines theorem applied to the LDE. The relation Eq. (2.11) 
shows that, to order-<5", the quantity 8r^ is directly obtained from the evaluation of Eren(0). As discussed in the 
introduction, we will use the Feynman rules described above to evaluate perturbatively the self-energy £ren(p) to 
order-5 4 from which we will get the nonperturbative values for k and TZ by using the PMS optimization procedure. 
The subscript "ren" in the self-energy means that this quantity also contains all diagrams which arise from the mass 
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counterterm vertex proportional to 5A$. For our purposes, the easiest way to obtain a perturbative expansion for 
(cA 2 )„ is to start from 



d 3 p N 



2 _|_ ^2 



(*>, 



p2 _|_ ^2 



(2.12) 



Like Sri , the order-<5" quantity (4> 2 )u^ is o btained by evaluating the self-energies to that order and subsequently 
expanding the series on the RHS of Eq. (2.12). Therefore, to obtain 1Z and k to order-<5 4 we need to consider the fifty 



four self-energy contributions shown in Fig. 1. 



III. EVALUATION OF r c TO 0{5 A ) 

According to the Hugenholtz-Pines theorem, Sri 4 ^ is obtained from the evaluation of all diagrams shown in Fig. 
1 with zero external momentum. To make this paper more pedagogical, let us do a step by step evaluation of r c 
up to order-<5 2 . To order-<5 one has only the tadpole contribution, a direct application of the Feynman rules for the 
interpolated theory and dimensional regularization (see appendix for more details) gives the finite contribution 



- SrP = EW (0) = -6u± + . (3.1) 
Carrying on to order S 2 one considers the contributions depicted by the first five diagrams of Fig. 1, which give 



<&f> = (0) = -Su 



Sn 



N ■ 



5 2 u 



V 
16tt 



N + 2 



-S 2 



16tt t) 



N ■ 



+ S 2 



1287T 



7V + 2 



-S 2 



u 2 (N + 2) 



(87T) 



18 



- +41n ( — ) 2.:-!<i4 

e \r] 



+ 5A S + 0(6 3 ) , 



(3.2) 



where M is an arbitrary MS mass scale 2 . Now, one replaces Sr c which appears at the right hand side with the 
value Sri^ obtained at the previous order so that the right hand side remains of orde r S 2 . Next, one sees that the 
setting sun, whose explicit evaluation follows those performed in the appendix [see Eq.( A18] )], displays an ultraviolet 
pole as e — * oo. In fact, within dimensional regularization, the only primitive ultraviolet divergence associated with 
the effective super-renormalizable three dimensional theory steams from the setting sun type of diagram with three 
internal propagators. The pole associated with this divergence fixes the mass counterterm coefficient in the modified 
minimal subtraction renormalization scheme, 



5 A s = S 2 



(N + 2)l 



18 



(3.3) 



As usual, this "vertex" must be considered also at higher orders (see Fig. 1) so diagrams whose divergences arise from 
"setting sun" sub-diagrams may be r endered finite pQ|. Now, it is easy to see how the "double scoop" contribution 
(fourth term on the RHS of Eq. ( |3.2| )) is exactly canceled due to the HP condition applied to r c at first order. One 
then gets the finite second order result 



- 6r^ = S( 2 ) (0) = -*„ JL 

07T 



N + 2 



5 2 u- 



rj (N + 2 



16tt 



u 2 (N + 2) 



(87T) 



18 



41n[M 

V 



2.394 



0(S 3 ) (3.4) 



Also at higher orders many contributions cancel. In special, any diagram with one or more tadpole sub-diagram(s), 
like the "double scoop" discussed at order-<5 2 , disappear. Then, the diagrams which really contribute to Sri^ are 



2 Note that our choice for the integral measure, Eq. ( [Al| ), has already taken care of constants like 4ir and e" IB which, otherwise, would 
appear in the setting sun logarithmic term. 
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those shown in Fig. 2, where one must consider the external lines as carrying zero momentum. At the same time, 
counterterm diagrams associated with the zero external momentum setting sun diagram (or sub-diagrams) could have 
been suppressed from that figure. However, we prefer to write them explicitly so that the same figure can be used 
again, facilitating the discussion in the next section. Using the results obtained in the appendix one obtains the fourth 
order result 



vr V 3 / 16tt V 3 / (4tt) 2 18 



In [ — | - o.r,!)77. r ) 

n 



, n fN + 2\ ,u 3 (N + 2) 2 r , r , u 3 (16 + 10A + A 2 ) r 

+ $ u-j- — ^— -6 — ) nn ,, 0-143848 + ^ —~ r, ,*, nn - 81.076 

64tt V 3 / r] 108(4tt) 3L 1 77 (4tt) 5 108 1 1 

,4 



(N + 2) (16 + ION + iV*) _ 5 ,uHN^f _ s ^N + *) [0A6m2] 

rf 6(4tt) 6 108 1 V 108(4tt) 3 1 1 18(4tt) 2 L J 

,4,(^ + 2), , .u 3 (l6 + 10A + A 2 ) r , ,u 3 (l6 + 10A + A 2 ) r 

iV .,. 0.0834 + S 4 — ± fA qiO.240 + S 4 — ± fA ^— ^ 30.31096 

18(4tt) 2L 1 t] (4tt) 5 108 1 J 77 (4tt) 5 108 1 1 



?7 2 (4tt) 6 648 1 J rf (4tt) 6 324 



The scale dependence of this quantity will be discussed in Section V. 



IV. EVALUATION OF (<f> 2 ) u TO ORDER S 4 

In principle, to obtain {(j) 2 )^ one should consider all contributions to the self-energy Erei(p) given by the diagrams 
of Fig. 1 with external momentum p. However, thanks to results of the previous section, one does not have to do the 
evaluation of all those graphs explicitly at this stage. In fact, one can immediately reduce the number of graphs to be 
considered by substituting the vertex Sr with the the appropriate critical value Sr c obtained in the previous section. 
Then, as for S^en(0), the set of diagrams which effectively contribute to Eren(p) reduces t o th ose shown i n Fi g. 2, 



but now one must consider the external lines as carrying momentum p. Substituting Eq. (3.5) into Eq. (2.12) one 

sees that the quantity which matters for the evaluation of {(j) 2 )^ is £rem(p) — Ercn(0). Diagramatically, this quantity 
is given by taking the graphs of Fig. 2 with zero external momentum and subtracting them from the same diagrams 
with external momentum p. This means that all diagrams which do not depend on the external momentum will not 
contribute in the evaluation of (tfi 2 ) at the critical point. For example, all the tadpole diagrams with any type of 
sub-diagrams will not contribute. As expected, the mass counterterm is a redundant quantity in the evaluation of 
{4> 2 )v^ because this quantity depends on the difference 

E^(p) - S (»)(0) = pW(p) + ^\p)} [Eft>(0) + E<; 1) (0)] , (4.1) 

where E^(p) is the divergent self-energy. For a general renormalizable theory, the quantity E^(p) represents all 
counterterms associated with the parameters of the theory (such as masses and coupling constants) as well as the 
wave-function counterterm associated with any eventual momentum dependent pole. At the same time, E^™ (0) 
involves the same counterterms except for the wave-function one. However, as we have already emphasized, in the 
three-dimensional case the only type of primitive divergence requires only a mass counterterm, which is the same for 
^div(p) an d Ejjj^(O). This means that in our case, E^(p) — E^(0) is always a finite quantity as shown explicitly in 
the appendix where it is also shown that this quantity is scale independent, as opposed to r c . Therefore, the type of 
diagrams which really matter for the evaluation of {(p 2 )^ are those shown in Fig. 3, which can be obtained expanding 



Eq. (2.12) to 0(S i ). Following the sequence of diagrams shown in Fig. 3 one can write 

d 3 p N f 5rf <5V *V <*V 



27r) 3 p 2 + r\ 2 \ p 2 + r\ 2 (p 2 + r/ 2 ) 2 (p 2 + ?? 2 ) 3 (p 2 + t] 2 ) 4 
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_ ^ [Ex^-IMO)] j3 2r; 2 [S 1 (p)-S 1 (0)] j3 [S 2 (p) - S 2 (Q)] 
p 2 + r/ 2 (p 2 + r/ 2 ) 2 p 2 + r) 2 

[S 3 (p) - E 3 (0)] _ J4 3^gi(g)-Si(0)] _ (54 2t ? 2 [S 2 (p) - E 2 (0)] 



p2 _|_ ^2 (p 2 +r] 2 ) 3 (p 2 + rj 2 ) 2 

§4 [S 4 (p)-S 4 (0)] _ J4 [S 7 (P)-S7(0)] _ g4 2r? 2 F3(p)-S 3 (0)] 

(p 2 + ry 2 ) (p 2 + ?7 2 ) (p 2 + ?7 2 ) 2 

J4 [Siq(p) - Sio(O)] _ j4 [S 5 (p)-S 5 (Q)] + ^ giCgj-Sig 

(p 2 + ?7 2 ) (p 2 + ?7 2 ) (p 2 + ?7 2 ) 2 

g4 [S 6 (p)-S 6 (Q)] _ j4 [S 8 (p)-S 8 (Q)] _ j4 [S 9 (p) - E 9 (0)] 



(p 2 + if) {p 2 + rf) (p 2 + rj 2 ) 

j4 P3 11 (p)-S„(0)] + 5 | 
(p 2 + ?T) J 

The details of the explicit evaluation of the Sj terms are given in the appendix. The final result we obtain is 

h 



2 Nrj t SNri t S 2 Nrj S 3 Nrj S 4 5Nri 



4tt 2 4tt 8 4tt 16 4tt 128 4tt 

^^#[0.143848] - ^^^[0.01168] - ^^[0.0610] 
7} 18(4tt) 3 l j 77 18(4tt) 3 l j 77 18(4tt) 3 l j 

^^ (16+1 .°: + ^ [8.0^] - ^S±|l [, 8270 „ 10 -3] 

^ ^ i™^, 30476] _ s ^mN+v mm 

._ +N 3 ) r , .u 4 N (44 + 32iV + 5iV 2 ) r 

^ ^ 3.15904767 - 5 4 ^r-—=± — qi.70959 
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V. NUMERICAL RESULTS FOR THE TEMPERATURE SHIFT 



In this section we will turn our, so far, perturbative evaluation into nonperturbative results using the PMS opti- 
mization prescription. Our analysis of results, including the selection of the relevant optima, will follow close ly t hose 
adopted in the applications which proved the convergence of this method for the anharmonic oscillator pi], [2^. g3|. ^4) . 
Some of the guidelines developed on those studies are essential for our present application. Let us start the optimiza- 
tion process with the scale independent quantity {4> 2 )u^ whose recent Monte Carlo estimate is {4> 2 ) u — —0.001198(17) u 

||. Before optimizing let us remark that all contributions to (t^ 2 )!™' are proportional to 5 n u n r] 1 ~ n and therefore, 
the PMS condition will imply solving a polynomial equation of degree n. As one may expect, many of those n roots 
which determine the optimum fj will be complex. Also, as observed in the anharmonic oscillator studies, most of the 
time the best results are in fact generated by the complex solutions |2l| . Since r\ is arbitrary we have no justification, 
a priori, to throw away its complex part. This means that our optimized physical quantities (<p 2 ) u and r c will have, 
eventually, complex parts whose meaning is to be interpreted according to the physics. Here, these two quantities 
are ultimately used to determine a strictly real physical quantity defined by the critical temperature. Therefore, for 
our purposes the complex parts of those two physical quantities are not relevant and will not be considered. Note 
that the imaginary parts of optimized physical observables have also been dropped in Ref. pl| where a different, but 
still valid, physical argument has been used. Finally, we shall follow the original PMS prescription and optimize 

{4' 2 )u^ and ri^ separately. This procedure was also adopted in the ultra relativistic case where it has produced good 
results p6| . 
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By truncating Eq. (4.3) to the first nontrivial order, order-<5 2 , setting 6 
two real roots 



1 and by applying the PMS, one gets the 



fj = ±0.0232332 u , (5.1) 

which give 

((j) 2 )W = =f0.002777326 u . (5.2) 

Applying the PMS to (<j) 2 ) { u 5] at order-<5 3 one obtains the following three solutions. The first, 77 = — 0.0475422 u gives 
((j> 2 )^ = 0.0045505 m while the other two 

fj = (0.0237711 ±0.0268995i) u , (5.3) 

yield 

(0 2 )(f) = -(0.00221912 ± 0.00150245i) u . (5.4) 

At order-5 4 one obtains the real solutions fj = 0.0439352 u which gives (4> 2 )^ — —0.00293974 u and fj = 
-0.0697993 u which gives ((f) 2 }^ = 0.00483554 m. The complex solutions are 

fj = (0.0129321 ± 0.04676942i) u , (5.5) 

from which one gets 

(0 2 )( 4 ) = -(0.00134323 ± 0.00213104i) u . (5.6) 

In order to select the appropriate roots we recur again to the AO convergence studies where the existence and 
behavior of optima families was fully investigated to order-J 47 [pl[ . There, it was observed that at a given order n 
each PMS solution belongs to a different family, the exception being complex conjugate solutions which belong to the 
same family. It was observed that, in the complex plane, the first member of a new family always lies on the real 
axis and also that a new family arises as n is increased by 2. Supposing that these findings may also be used in our 
three-dimensional problem, we may identify two families whose first members lie on the real axis at order-<5 2 . Family 
1 starts with the positive real solution fj = 0.0232332 u and family 2 with the negative real solution fj = —0.0232332 u. 
No new families arise when one goes to the next order and the real negative solution fj = —0.0475422 m is just another 
member of the family of negative real solutions (2) while the complex conjugate optima with positive real parts 
fj = (0.0237711 ± 0.0268995?) u are taken as belonging to family 1. At order-<5 4 , family 2 gets another member given 
by fj = -0.06983 m, whereas family 1 gets fj = (0.0129321 ± 0.04676942i) u. 

As we have increased the order by 2 one effectively sees the appearance of a new family whose first member lies 
on the real axis and is given by fj = 0.0439352 u. We can now roughly examine the convergence of our results. The 
values obtained with the optima belonging to family 1 are {4> 2 )u ] = -0.002777326 u, Re[(</> 2 )1 3) ] = -0.00221912m 
and Re[(< ? !) 2 )l 4) ] = -0.00134323 u. Family 2 gives {(j) 2 )^ = 0.002777326 u, (0 2 )i 3) = 0.00405505 u and (0 2 )1 4) = 

0.00483554 u, whereas family 3 gives (</> 2 )! 4) = -0.00293974 u. Note that the first (0 2 )„ value predicted by family 
3 is only about 5% greater than the first value predicted by family 1. It is very likely that family 3 will become 
complex and, as for the AO, as we go to higher orders families 1 and 3 will predict very similar values converging to 
the exact value. Family 2, on the other hand, seems to have only real components. It predicts values of (4> 2 ) u which 
increase order by order with a sign which is opposite to the one predicted by families 1 and 3. Moreover, in the AO, it 
was observed that the complex families have better convergence behavior than the purely real families. This analogy 
indicates that family 1 should produce converging results. 

We can justify pushing the analogy in between our effective three dimensional model and its one dimensional version 
that far by remarking that, at least to the order we consider here, {4> 2 )u^ can be expressed as a power expansion of 
the form 
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(^ 2 >i 4) = at - *) 1/2 r^ , (5.7) 

i=0 

where B a ~ 10 _1 , Si = 0, B 2 ~ 10~ 5 , S 3 - 1(T 6 and B 4 - 10~ 7 . This structure is similar to the one found in the 
one dimensional case. This hints that both models may have similar convergence properties making our procedure 
more legitimate. It is also worth pointing out that in our previous work, Ref. Q, we had only the order-5 2 result 
and it was not possible to do the same type of comparison among the solutions to find an acceptable pattern of order 
by order corrections. There, to choose among the two possible solutions, (0 2 )L 2 ^ = T0-002777326 u , we had to use 
different arguments and were also guided by results found with other methods. By considering higher orders, as we 
have done here, we can overcome this problem and the negative result, (4> 2 )^ = —0.002777326 u , naturally appears 
as the one which belongs to the most well behaved sequence of order by order corrections. 



We are now in position to evaluate un — A^ 2 )^, so that c\ can be determined via Eq. ( |1.7[ ) with the optima 



contained in family 1. As one could expect, fj is always proportional to u since the latter quantity is the only 

quantity with mass dimensions appearing in {<j) 2 )u^ ■ This means that the optimum value for the non-interacting 

vacuum expectation value (<fi 2 )^ will be zero at any order. This agrees with the results of Ref. ||, where it was 
shown that this is indeed the value obtained when the theory is regularized with dimensional regularization. Then, 

uk, = A(0 2 )c™' ) = (0 2 )!"^ from which one finally obtains c\ = 3.06, c\ — 2A5 and c\ — 1.48 at order 5 2 , 6 3 and 8 i , 
respectively 3 . These results are compared with other analytical and numerical results in Table 1. 

It is instructive to examine the topology of the diagrams contributing at each order so that we can establish the links 
with other nonperturbative methods. At second order the non-trivial contribution arises from the setting sun (one 
plain bubble) type of diagram. At third order one has, besides the setting suns with insertions, a new contribution 
which arises from the two plain bubble type of diagram (ninth graph shown in Fig. 3). However, this contribution, 
belongs with the setting sun to a class of diagrams that would appear in a plain bubble sum or in the leading order of 
a 1/N type of calculation. At fourth order one considers again a three plain bubble contribution (eighteenth diagram 
of Fig. 3) but more radical changes arise via other type of vertex corrections like the correction to the plain bubble 
that comes from the nineteenth and twentieth diagrams of Fig. 3. Finally, the last diagram contains a different type 
of vertex correction that would appear on a ladder type of summation. In fact, one can easily evaluate which are the 
individual contributions of the five-loop diagrams shown in Fig. 3. The first of them gives a contribution (in unities 
of u 4 /rf) of approximately 1.9 x 10 -9 , the second gives 2.8 x 10~ 8 , the third 2.7 x 10 -8 , the fourth 8.86 x 10 -9 and 
the fifth gives 2.6 x 10 -8 . These numbers show that, at this order, the total contribution from the ladder and bubble 
correction type of contributions (third and fifth) are effectively twice that of the plain three bubble one. 

It is also easy to see by drawing that the only corrections which may appear at odd orders are those due to the 
doubling of a bubble that already appeared at the previous order (increasing the "bubble chain"). At the same time, 
at even orders, one is allowed to insert a new bubble anywhere creating diagrams with completely different topologies. 
In other words, in a perturbative expansion of {4> 2 )u \ new topological classes of graphs can arise only at even orders. 

One can now appreciate that the reason our order-<5 2 result c% = 3.06 Q, obtained by optimizing only one setting 
sun contribution, compares so well with the value c\ = 2.90, found by resumming setting sun contributions in a 
self-consistent way 0, is a consequence that both approximations consider the same type of diagrams. On the other 
hand, when going to order-5 3 one considers a new diagram but which, together with the setting sun, would also be 
considered in a large-iV calculation. In the LED, its effect is to reduce the second order result to c% = 2.45. Let 
us consider, for the moment, the only order-5 4 contribution that would also be considered in a large-iV calculation. 
Graphically this contribution is displayed by the second of the five-loop terms in Fig. 3. Not surprisingly, we obtain 
the value c\ = 2.32 which is very close to the ci = 2.33 value obtained with the 1/N method at leading order g and 
the numerical differences may be due to the fact that we have considered our symmetry factors in full, not only the 
highest power of N. The four remaining five-loop contributions would be considered in an 1/iV type of calculation 
to the next order. Such a calculation has been performed by Arnold and Tomasik [^| who found c\ = 1.71, which is 
approximately 27% smaller than the leading order result. In our case this fact is confirmed at order-5 4 , where the net 
effect of considering diagrams which would belong to a next to leading order 1 /N evaluation is to decrease the value 
c\ = 2.32 obtained with the graph that would appear at leading order in the same approximation by roughly 35%. As 
before the numerical differences must be due to the full consideration of powers of N in each symmetry factor. It is 



3 At this stage it should be clear that it is preferable to optimize (<f> 2 )u^ rather than A^ 2 )^' because the latter quantity is less r) 
dependent. 
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not our aim to establish here a formal relationship among the different approximations. Nevertheless, the discussion 
above can serve as a guide to understand how the LDE captures part of the nonperturbative physics contained within 
the SCR and 1/N approximations. 

In order to evaluate the coefficient c' 2 ' we now turn to the optimization of the scale dependent r c . Setting 5=1 
and applying the PMS to r c generates one positive, real optimum given by fj = u/6tt. It is important to note that 
this PMS solution is a scale independent quantity. In fact, r c depends on the (MS) mass scale through the term 
proportional to u 2 ln(Af /rj) which appears in the order-(5 2 setting sun term. It is then easy to see that when this term 
is derived with respect to r\ the scale dependence automatically disappears turning our optimization procedure into a 
scale independent process. As discussed below this situation will be verified at any order in 5. 

Next, in order to get a numerical result for the optimized r c one must fix a scale and here we choose M = u/3 
which is the same scale 4 used by Arnold, Moore and Tomasik in Ref. ||] where the result found for this quantity 
is r c (M = u/3) = 0.001920(2) u 2 . The relatio n in between the values of ri S \ evaluated at two different (MS) mass 



scales Mi and M2, can be obtained from Eq. (3.5) and reads 



^ 4) (Mi) rj 4 \M 2 ) (N + 2) (MA 

It is not too difficult to see that this relation will be verified at any order in 6. At order-5 2 the only diagram which 
is scale dependent is the setting sun. At a higher order (n > 3) this order-<5 2 contribution can only appear as a 
subdiagram. At the same order a similar graph appears, but this time 5r replaces the setting sun insertion. However, 
the "vertex" Sr c is always replaced (see Sec. Ill) by its expansion in S which contains, at order-<5 2 , exactly the same 
scale dependent term as given by the setting sun, with a reversed sign. This means that, apart from the order-5 2 
setting sun, all contributions to Sr^ are automatically scale independent. Optimizing our order-<5 2 result one gets 

r (2) = 0.00315236 u 2 . (5.9) 
Proceeding to next order the PMS gives two complex solutions, fj = (0.0353678 ± 0.0550091i)u 2 which yield 



r 



(3) _ 



(0.00221321 ± 0.0000966H) u 2 ■ (5.10) 



Finally, the order-5 4 optimization results are the real solution 77 = 0.0659334 m, which yields = 0.00246153 u 2 and 
the complex solutions fj = (0.00947463 ± 0.0797262i)it 2 , which generates 



r< 4) = (0.00165411 ±0.000772567i)u 2 . (5.11) 

As in the previous case one sees that the first optima family starts with a real value at order-<5 2 and turns into 
a complex family at order-5 3 . At order-J 4 it receives a new complex member. The first family generates the real 
values r { c 2) = 0.00315236m 2 , Re[r£ 3) ] = 0.00221321 u 2 and Re[r< 4) ] = 0.00165411 u 2 which are our selected values. 
Then, using Eq ([01) together with the optima values obtained for k and 1Z we obtain, order after order, the results 
c'2 = 101.4, c'2 = 98.2 and c' 2 ' = 82.9 for the order-a 2 nonperturbative coefficient. As for ci, these results compare well 
with the Monte Carlo estimate, c' 2 ' = 75.7 ± 0.4. 



VI. CONCLUSIONS 



We have used the linear 6 expansion to evaluate nonperturbatively the numerical coefficients appearing in the 
expansion for the transition temperature for a dilute, homogeneous, three-dimensional Bose gas given by T c = Tq{1 + 
Cian 1 / 3 + [c' 2 h^an 1 / 3 ) + C2]a 2 n 2 / 3 + 0(a 3 n)}, where T is the result for an ideal gas, a is the s-wave scattering length 
and n is the number density. This expansion for T c incorporates the effects of non-zero Matsubara modes |8|, [| . While 
the coefficient c' 2 has been exactly evaluated using perturbation theory the question about the numerical values of the 



4 Our notation for the mass scale (M) is different from the one used by the authors in Ref. JfJ (M). 



13 



other two coefficients, c\ and c' 2 ' remains open and has been the object of recent investigations. The reason behind 
this difficulty is the fact that these coefficients can only be obtained in a nonperturbative way. 

Due to the Hugenholtz-Pines theorem the first non-trivial contribution appears at an order where one has to 
consider, at least, momentum dependent two-loop self-energy diagrams. Considering higher order terms, so as to get 
more accurate results, becomes rapidly difficult within the existing nonperturbative methods as discussed in Ref. [[1, 
where the authors state that the complexity of the mathematical problem does not allow a definitive prediction of 
the prefactor c\, of the term linear in a, from an analytic analysis. On the other hand, two recent numerical results 
obtained with lattice simulations, which predict c\ ~ 1.30 H, |§, are being taken very seriously. In a previous 
work, Ref. Q, we have applied the LDE to this problem obtaining the value c\ ~ 3.06 at the first non-trivial order 
(<5 2 ). However, the quality of that application was difficult to infer, from a quantitative point of view, since only one 
approximant had been used. On the hand, the fact that at order-5 2 with only one graph the optimization procedure 
was able to generate a result numerically similar to the one obtained with a self-consistent resummation (SCR) of 
two-loop momentum dependent contributions was encouraging. At that time, we were not in position to elaborate 
any further about the convergence behavior of that result. 

In the present work, we have again explicitly shown that the LDE method offers, as its major advantage, the 
possibility to select, evaluate and renormalize a physical quantity exactly as in the familiar perturbative framework. 
Here, the contributions appearing at each order are not selected according to their topology as within most nonper- 
turbative analytical cases. Contrary to some previous unfounded criticisms, no uncontrolled errors arise in this type 
of perturbative calculation, most notably in this application, where even the most cumbersome five-loop contribu- 
tions have been fully considered and evaluated without recurring to any approximations as shown in the appendix. 
Another advantage is that one does not have to worry about infrared divergences, since, during the formal evaluation 
of graphs, the LDE arbitrary parameter naturally acts as such before disappearing during the optimization process. 
Also, the fact that a convergence proof for the quantum mechanical analogue of the model considered here does exist 
pl| . p2| , p3| , pi[ is an extra bonus. 

At first one could think that the multiplicity of possible real and complex results generated by the PMS constitutes 
the most serious disadvantage of the LDE. Nevertheless, the quantum mechanical convergence studies of Ref. [El), 
have shown how meaningful nonperturbative physical results can still be obtained. As discussed in the text, those 
studies have been crucial to our application for some important reasons like showing how the possible solutions gather 
into real and complex families and emphasizing that better results are generated by the complex ones. We recall that, 
although different physical arguments have been used in each case, the imaginary parts of the optimized physical 
observables generated by the complex families have also been dropped out in Ref. pl|. As already mentioned, 

our effective model displays the same series structure for the physical observable (4> 2 }u cis its quantum mechanical 
counterpart. Taking all these facts into account we were able to obtain the results c\ = 3.06, c\ — 2A7 and c\ = 1.48 
at second, third and fourth orders, respectively. Our results approach, order after order, the recent Monte Carlo 
estimate, c\ ~ 1.3. 

Comparing our results and the topology of the diagrams considered here with those belonging to the self-consistent 
resummation of setting suns (SCR) and the 1/N approximation at leading (1/A-LO) and next to leading (1/iV-NLO) 
orders we made clear that our results are not a mere coincidence. In fact, the PMS is successively chopping, order 
after order, nonperturbative information contained in those approximations. Our results confirm the decrease in the 
value of Ci observed successively with the SCR, 1/N-LO and 1/A-NLO. The numerical differences may be due to the 
fact that we do not make any distinction among the different powers of N which appear on the symmetry factors 
since the LDE was envisaged to cope with arbitrary N. 

We remark that a problem regarding the sign of the coefficient c\ , which appeared in our previous application, has 
disappeared at this higher order evaluation. We have also investigated the quantity r c by evaluating all self-energy 
contributions, with zero external momentum, up to order-<5 . Once this quantity was optimized we have obtained the 
values c' 2 ' = 101.4, c'2 = 98.2 and c'2 = 82.9 for the next nonperturbative coefficient at second, third and fourth orders, 
respectively. These results are in good numerical agreement with the Monte Carlo result, c 2 ' — 75.7 ||. 

In summary our analytical investigation seems to support, order by order, the results obtained with other three 
analytical nonperturbative methods. Our fourth-order numerical results compare well with the recent results found 
in Refs. ^|, |). Additionally, there is an exciting possibility that the method may offer a way of making a definitive 
analytical prediction for the nonperturbative coefficients c\ and c'2, which we are currently investigating. 
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APPENDIX A: EVALUATING THE HIGHER LOOP TERMS 

To make this work self-contained we shall outline, in this appendix, the details of the explicit evaluation of all 
Feynman diagrams considered in the evaluation of (4* 2 )^ for arbitrary N. We also remark that working out symmetry 
factors for many loop contributions with generic N is a problem on its own. Here, we have used the methods developed 
by Kleinert's group in Berlin [^7| . 

We regularize all diagrams with dimensional regularization in arbitrary dimensions d = 3 — 2e and carry the 
renormalization with the modified minimal subtraction scheme (MS). So the momentum integrals are replaced by 
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where M is an arbitrary mass scale and — 0.5772 is the Euler-Mascheroni constant. Very often, in evaluating the 
contributions to (</> 2 ) u one considers the integral 
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(A2) 



1 



pi _|_ 



47T 



1 + 

ze! 

4 



2 In 



2 In 



M 
V 

M 
2rj 



2 - ln(4) 



41n 



+ 0(e 3 ) 



and its derivatives with respect to i] 2 



1 



1 



(p 2 + rj 2 ) n (n — 1)! \ dr\ 



1 



p2 _|_ ^2 



(A3) 



(A4) 



Let us now consider the three-loop contributions to {4> 2 ) u with any number of, external and/or internal, Srj 2 
insertions. Their general form is 



(<V) ^ 



(p 2 + tfy 
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where c is defined below and n determines the number of external (to the setting sun) Srj 2 insertions. At the same 
time the insertions, internal to the setting sun, are taken into account by 
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where M. defines the multiplicity of equivalent internal 5r] 2 insertions. This general contribution to ((j> 2 ) u can be 



written as 



(5n 2 ) n - J N 



5 c zZ a (p) _ xn+m+ j + h-3 N(N + 2)M 2 . 2 \n+m+j+h-5 
(p2 +n 2)n ~ 18 W J 
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where c — m + j + h— 1 labels the order of the two-loop (setting sun) self-energy term. Now, we can merge all 
propagators through the use of standard Feynman parametrization, given as usual by 



1 _T[x + y] f 1 ^ a x - 1 (l-a) y - 1 
a x by ~ T[x]T[y] J ° [aa + 6(1 - a)] x +y 



d^i — rrn W> (*,»><>). (A8) 



or other generalizations. One then gets 

5 c Y> a {p) 5 n+m+j+h-3 N ( N + 2 )M u 2 T[n + m + j + h- 9/2 + 3e] 



(5r ] 2 y i - 2 N 



p (p 2 + V 2 ) n (4tt) 9 / 2 18 r] T[n}T[m}r[j}T[h] 

e^M 2 \ 3e [\ , a , g(a)g(P)g{i) 



where 



1? J J dad ^ [g{ ^ M)] n +m+]+h -W > ( A9 ) 

g(a) = qJ'-^I - af-^ail - a )]-J-h+3/2-c , (A10) 

g(j3) = /^(l - l3y +h - 5 / 2+e [l3(l - p)]-j-h-m+3-2e ^ (An) 



and 



g(a ;/ 3 ;7 )^ 7 +^+ 1 7 . (A13) 
1 — p po;(l — a) 

Then, for given n, m, j and ft one performs the expansion in e keeping the poles and finite terms as usual. For most 
situations found in the present work the integrals over the Feynman parameters need to be evaluated numerically. 
Here we use Monte Carlo and Vegas techniques to perform those integrations. We have taken particular care to keep 
the numerical errors less than approximately 1% in our final numerical results. 

One must be careful in carrying out the e expansion in the expression above since sometimes the divergences can be 
hidden on the exponents of the Feynman parameters. Since m, n, j and h are positive integers (n > 2, m, j, h > 1) one 
sees that 5(7) has a pole as e — > when j = h = m = 1 corresponding to a setting sun diagram without internal Srj 2 
insertions. This is the only situation where one has ultraviolet divergences for these contributions. For j = h = m = 1 
the actual divergence appears in the term (1 — 7) 2c ~ 1 contained in 5(7) and it will appear asal/e pole if one integrates 



/' 

Jo 



1 dl (A14) 



by parts. The case p = follows essentially the same lines and the general result for the setting sun type of contribution 
with any internal and/or external Sr] 2 insertions is 

2 X a-2 N f ^q(O) 5 n + m ^+ h ~ 3 N(N + 2)M 2 r[m + j + h-3 + 2e] 

[ V> J P (p 2 +V 2 ) n ~ (4tt) 3 18 ^' T[m}r[j}T[h] 

?7 2 / Jo M [ff(a,/?)]"*+-H- h - 3 + 2e J p (p 2 + r 1 2 y K ' 

where the integral over p can be readily obtained by one of the methods discussed in the beginning of this appendix 
and 



16 



g(a,P) = 



1-/3 0a(l-a) 



(A16) 



Note that the gamma function in Eq. (A15) displays an ultraviolet pole when m = j = h = 1 and is finite otherwise. 

The first contribution (sixth diagram of Fig. 3) of this type appears at 0(S 2 ) with n = 2, m = j = h = 1. Since 
there is just one graph like this, A4 = 1, and one writes 



N 



(5 2 Si(p) ^N(N + 2)u 2 
= o 



l P (p 2 + v 2 ) 2 

The p = contribution is given by 



18(4tt) 



T) 



1 , M 
- + 6 In — 

e V. 



- 4.93147 



(A17) 



JV 



^ 2 Si(0) 

p {p 2 + v 2 ) 2 



18(4tt) 3 r] 



6 In 



3.78069 



(A18) 



The l ast t wo equatio ns re produce the results found analytically by Braaten and Nieto in Ref. |38| . Note that although 
Eq. (A17) and Eq. (A18) diverge, their sum is finite and scale independent. Together, they give the contribution 



N 



3 2 [Si(p)-Si(0)] 
{p 2 + r\ 2 ) 2 



,u 2 N(N + 2) 
V 18(4tt) 3 



0.143848] 



(A19) 



which is exactly the result found in our previous work, Ref. |4|. Now, we turn to the evaluation the setting suns with 
insertions. The most expedient way would be to do the replacement rj — > 77(1 — 8) 1 / 2 and then expand the squared 
root to the desired order. However, this procedure would only give the total contribution at each order. In order to 
have absolute control about each single contribution we prefer to use our general expressions Eqs. ( |A9D and ( A15| ). 
We have checked both procedures finding that they agree to each other within 1% which is reassuring since when 
obtaining the diagrams with insertions via 77 — > 77(1 — 8) 1 / 2 one has a result which may be considered exact since this 
expansion starts from Eq. ( Al9| ), and this result agre es with the analytical results of Refs. ||, [58|. Also, at order-<5 4 , 
the general relations given by Eqs. (A.9) and ( AL5 ) have proven to be very useful in the evaluation of diagrams 
containing the setting sun as insertion. For the eighth graph of figure 3 one has three cases similar to n = m = 2 and 
j = h=l (M = 3). This gives 



N 



S 3 E 2 (p) 



(p 2 + 7] 2 ) 2 



= 6" 



u 2 N(N + : 
V 18(4tt) 3 



-[0.188] , 



(A20) 



and 



N 



£ 3 E 2 (0) 
(p 2 + rj 2 ) 2 



7] 18(4tt) 3 1 1 



(A21) 



which lead to 



A / = S^ N [ N + ? [0.0610] . 



(p 2 + i 1 2 ) 2 ry 18(4^) 3 

For the seventh graph of Fig. 3, n = 3 and m = h = j = 1 with M = 1, one has 



(A22) 



N 



<5V£i(jj) ,u 2 N(N + 2) 1 



P (p 2 + V 2 ) 3 



= 5 



rj 18(4tt) 3 32 



1 M , 

- +6 In — I - 1.968 
e \r/ 



(A23) 



and 
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N 



<5V£i(0) _ x3 u 2 N(N + 2) 1 



p (p 2 + iff t] 18(4tt) 3 32 



6 In 



M 



1.781 



(A24) 



which lead to 



v /ggfeM^M = -^m+g [0 . 01168] 



(p 2 + ?7 2 ) 3 



7? 18(4tt) 3 



(A25) 



where the factor of 2 on the RHS accounts for the two possibilities of external insertions (see Eq. (4.2). For n = 4 
and m = h = j = 1 (tenth diagram of Fig. 3) one gets, with JA = 1, 



-N 



<jVSi(p) _ j4 iv(iv + 2) i m 2 



(p2 +r? 2)4 18(4tt) 3 64 ry 



-+61n( — ) 1.2011 1 

e V 7? 



(A26) 



and 



which leads to 



. f gVgi(0) _ 4 JV(JV + 2) 1 77 2 

7o (P 2 + »? 2 ) 4 ~ 18(4tt) 3 64 7? 



i + 61n( — ) - f/l -KKs 
e V V 



(A27) 



A / 3*Yfi(p)-Si(0)] = _,4gg+g)^ [2 . 8270 x 10 - 3] 



(p 2 + 7? 2 ) 4 



18(4tt) 3 77 



(A28) 



where the factor of 3 on the RHS accounts for the possibilities of external insertions (see Eq. (4.2). There are three 



cases (A4 = 3) similar to the case 77 = 3, m — 2 and j = h = 1 displayed by the eleventh graph of Fig. 3. One gets 



N 



p (p 2 + V 2 ) 3 18(4tt) 3 77 



[0.0586] , 



(A29) 



and 



, <S 4 77 2 £ 2 (0) ,N(N + 2) u 2 r 

N I -4 — = S 4 A,. J — [0.0650] 



(p 2 +7/ 2 ) 3 



18(4tt) 3 tj 



(A30) 



which lead to 



N 



5V[2S a (p)-53 a (0)] s4 N(N + 2)u 2 



(p 2 + 77 2 ) 3 



18(4tt) 3 77 



[7.7318 x 10- 



(A31) 



The twelveth graph of Fig. 3 has .M=3,T7 = 2,m = 3 and j = h = 1 leading to 



N 



J^ =j «™> ^[0.058638] 
{p 2 + rj 2 ) 2 18(4tt) 3 77 L J 



(A32) 



and 



N 



5 4 S 4 (0) = __ 54 A^(A^ + 2)7i 2 



„ (p 2 + 77 2 ) 2 18(4^) 3 77 



[0.083246] , 



(A33) 



which give 
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Finally, for the thirteenth graph, M = 3, n = m = j = 2 and h = 1 from which one gets 



^(P) -^^+^[0.0234552], (A35) 



( p2+r? 2)2 18(4tt) 3 77 



and 



N j Jp% . _ 4 <ffi±a>Hi |( ,0417] , (A36) 
' (p 2 + ?7 2 ) 2 18(4tt) 4 r? 



which gives 



Let us now consider a general four-loop contribution with any number of internal and/or external 5r] 2 insertions. 
After performing few shifts on the integration variables one gets 

J P ( P 2 + V 2 ) n ' 108 U ' 



mk t (t 2 + v 2 ) n (q 2 + v 2 ) m (k 2 + v 2 Yl(p + q) 2 + v 2 ] h 

' (A38) 



[(p + fc) 2 + f] 2 Y[(p + 1) 2 + v 2 V ' 



where d = m + l + h + i+j — 2 labels the order of the three-loop self-energy term. Then, proceeding as in the three 
loop case one finds 



. „ _ 2 ■ 5 d Z b (p) _ Zn+m+l+h+i+j-A u 3 M ( lg + WN + jy2) f M 2 exp ^ 

[Orj ) iv / — 



' p (p 2 +?7 2 )™ (4tt) 6 ?7 2 108 V V 2 

T(n + m + I + h + i + j - 6 + 4e) 



4 c 



X 



T(n)T(m)T(l)T(h)T(i)T(j) 



x / *^W r ,, 6+4 , (A39) 



o 



where 



/(a) = a 1 ' 2 - 1 -^! - a) 1 / 2 -'- £ , (A40) 

/(/?) - ^-^(l - /3) 1 / 2 - m - £ , (A41) 

/(7) = 7 l/2- J -- e(1 _ 7) l/2-„-e 5 (A42) 

/(6>) = r i+ ' l - 5 / 2+e (l - 6»)"+^- 5 / 2+e , (A43) 
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f((j)) = _ ^ 

and 

/(Q , A7 , M ,= * + ^ + (lziMW). (A45) 

a(l-a) /3(l-/3) 7(1-7) 

As far as renormalization is concerned one should note that those type of four-loop contributions are always finite. 
The four-loop contribution whose self-energy has zero external momentum reads 

where 

E 6 (0) = NS 3 u 3 M ^ + ™ N + ^ (^ 2 )»- 2 (Vr +Wl+t+J - 5 

108 



1 



21m 



gfc* (<? 2 + ?? 2 ) J (fc 2 + ?7 2 )'(i 2 + v 2 ) h i(q + k) 2 + v 2 Yi(q + 1) 2 + V 

Proceeding as above one gets 

S d M0) xn+m+l+h+i+ , 4 „u 3 (v 2 ) n - 5/2 M (16 + 10JV + AT 



(A47) 



(S^Y^N / / 9 Sb( °j = S n+m+l+h+i+j-i N l 

Jv vp +v r 



2\ 



lp{p 2 + V 2 ) n (4tt) 9 / 2 108 

' M 2 exp7 B \ 3e r(m + I + h + i + j - 9/2 + 3e) 



r? 2 J r(m)r(z)r(ft)r(i)r (7 i 

x /"cmw* 7/(«)/W(7.«) 



o 



X / / 9 1 9N„ » ( A48 ) 



where /(a) and are given by Eqs. ( A40 ) and ( A41 ). Also, one has 



/( 7 ,0) = (1 - 7)^M7(1 - e)] l+l+c - 5 ' 2 K1 8) h+m - 5 ' 2+t , (A49) 

and 

/(a, (3, 7, 0) = (1 - 7 ) + + -577^ ■ (A5Q) 
a(l — a) p(l — p) 

The first four-loop contribution of this type appears at order-5 3 and is displayed by the ninth graph of Fig. 3 and it 
has n = 2, m = l = h = i= j = l and A4 = 1. The contributions to {<p 2 ) u are given by 

- Pn f T Il% . -^j y i6+1 °" + W2 } [32 .4388] . (A51) 

7 p (p 2 + ?7 2 ) 2 ?7 2 (4tt) 6 108 1 j y ' 

and 

^ / . *g> . ,.4 » £tt"£±^ [40 . 538] , (A52) 

Jp(p 2 + V 2 ) 2 V 2 (4tt) 6 108 L J v ' 
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which lead to 

Jp (p 2 + V 2 ) 2 V 2 (4tt) 6 108 1 J V ; 

As for the three-loop case one could use the equation above to obtain a scries expansion which would give the total 
contribution of graphs with insertions to any order in 5. However, we prefer to perform the individual evaluation 
of each contribution in order to achieve more control over the series expansion. At order-<5 4 the first contribution is 
displayed by the fourteenth graph with M. — 1, n — 3 and m = l = h = i= j = l. One then obtains 

and 

s' N [ = jl ^±m±n mm , (A55) 



p ( P 2 + V 2 ) 3 ii 2 {Anf 108 



which lead to 



~ J P (P 2 + V 2 ) 3 " ^ 1O8 [ °- 85984] ' ( A56) 

where, once more, the factor of 2 accounts for the two possibilities of internal insertions in accordance with Eq. (|f .2]). 

Next, let us consider the case illustrated by the fifteenth graph of Fig. 3, which has M. = 1, n = j = 2 and 
h = i = l = m=l. After evaluating the integrals one gets 

_ rtv / T ^4- . » (" i+1 | " + JV2 > [3,8221] , ,A57) 

7 P (p 2 + ?7 2 ) 2 v (47r) 6 108 



and 



^ f . 4JL <' c + » + "') [, 12 ,| , 

7 P + ?7 ) (4tt) 6 108 L J 



(A58) 



which lead to 



_ S , N j M-yi _ ^(»tMt»i w . (A59) 

J P ip + V ) V (4tt) 6 108 

The remaining four-loop contributions to this order are evaluated using the case displayed by the sixteenth diagram 
of Fig. 3 which has M. — A, n = m = 2, l = h = i = j = l and whose result is given by 

f <5 4 E 5 («) W N (16 + 107V + iV 2 ) , 

N / -, = byP ' = -^-y—jj i — ^ 9.85072 , A60 

J P {p +V ) V ( 47r ) 108 



and 



w / . ,4 jl (*±m±^ [15 , 5548] , (A61) 



which lead to 
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7 P (p 2 +7? 2 ) 2 ?7 2 (4tt) 6 108 1 J 

Let us now consider the five-loop contributions. The first one is given by the seventeenth graph of Fig. 3, 

fi„ ( P'W-y . [0.87339] , 

7 P (p 2 + ?? 2 ) 3 ?7 3 (18) 2 (4tt) 7 

where the individual contributions are given by three terms starting with 



5 A N 



I 



Fi(p)] 2 

(p 2 + r] 2 ) 3 



4 u 4 N(N + 2) 2 r(3/2 + 5e) / M 2 cx P7g y t 
ry 3 (18) 2 r(3)(47r) 15 /2 ^ ^2 J 

h(a)h{0)h(y)h(0)h{<f>)h( X ) 



x / dad/3djd9d(j)dx 



1 



[Ma,/3,7,M,x)] 3/2+5e 
3 <f>h(a,f3,j,e, X ) 9 07[l-ft(a,)9)] 



4h(a,f3,f,6,(f>,x) 8 /i(a, /3, 7, 0, </>, x) 



15 2 7[1 - ft(a, /?)]/i(a, /?, 7, 0, x) 
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+ 



+ 



[/i(a,/3,7,0,<?!>,x)] 2 
5 #(a,/3, 7 ,0,x) 18 ^[1-/1(0;,/?)] 



2 ft(a, /3, 7, 0, 0, x) 4 ft(a, /3, 7, 0, 4>, x) 

8 [/i(a,/3, 7 ,^0,x)] 2 
5 07[l-/i(a,/3)] 25 2 7 [1- ft(a,/3)]ft(a,/3,7,0,x) 



2/i(a,/3, 7 ,M,x) 



where 



fc(a) = [a(l - a)]" 1 / 2 " 6 , 
M/3) = (l-/3)" 1/2+e [/3(l-/3)]" 2 
M7)-7(l-7) 2e , 

/>(0) = [6(1 - 0)]-V2-« , 

M0)-0 1+2e (i-^>) 2e , 
Mx) = (i-x)- 1/2+£ [x(i-x)]- 2 

/3) = — !-r + -5—^ r , 

1 — p pa(l — a) 



(A64) 

(A65) 
(A66) 
(A67) 
(A68) 
(A69) 
(A70) 

(A71) 



22 



h(fi, X ) = 



1 



1 



i-x xO(i-o) 



(A72) 



fc(a, /3, 7, 0, X ) = 7 - h(9, x) + (1 - l)h(a, 0) 



and 



(A73) 



h(a, 0, 7 , 0, <l>, X ) = l<t> + (1 - 0)M*> *) + 0(1 - 7)M«, /3) 
After performing the expansion in e and integrating numerically one obtains 



(A74) 



S 4 N 



Fi(p)] 2 

P (p 2 +?7 2 ) 3 



f4 u 4 iV(iV + 2) 2 



?7 3 1296(8tt) 



10 In 



71/ 



1 



4.419 



301n z 



,1/ 



(44.19) In 



M 



Now, expanding Eq. ( A9), with n = 3,m = j = h = l to order e and considering 



20.0158 



(A75) 



Ei(0) = -S 2 



one gets 



u 2 (N + 2) f 1 
8^ 18 - 



41n 



A/ 



2.3911 + e 



;ln 2 



r? 7 3 



9.5644 In 



M 



4.3127 



, (A76) 



p (p 2 +?7 2 ) 3 



-5 



4 w 4 /V(/V + 2) 2 f 1 



?7 3 648(8tt) 5 i <•- 



10 In 



71/ 



4.3067 



50 In 2 f — ] - (42.937) In f — ) + 18.9756 }> .(A77) 



The final contribution to this d iagram is obtained by considering Eq. ( |A2j ), with n — 3, expanded to order e 2 and by 
taking the square of Eq. ( A76 ) which leads to 



S 4 N 



Pi(0)] 2 _ x4 u 4 N(N + 2) 2 f 1 



p (p2 + 7? 2)3 ^3 i296(8tt) 



5 1 ,2 



10 In | — ) -4.1684 
V J 



50 In" ( - (41.684) In 



18.6434 



Next, let us consider the eighteenth graph of Fig. 3 whose contribution comes from 



S 4 N 



S 6 (P) 
(p 2 + 7y 2 ) 2 



= 6 4 



N (40 + 32iV + 8N 2 + N 3 ) 



r(3/2 + 5e) 



rf (4tt) 15 / 2 648 

1 dadf3d 7 d9d<t>d X d( ^M/%(7M%(^(xM0 



M 2 exp7s 

9 



5e 







[y(a,/?, 7 ,M,X,C)] 3/2+5e 



where 



(A78) 



(A79) 



y(a) = a[a(l - a)]" 3 / 2 " 6 , 



(A80) 



y((3) = [0(1 - )9)]-Va-« , 



(A81) 



y( 7 ) = [ 7 (l- 7 )]-V2-e ) 



y(9) = [6(1 - O)}- 1 ^ , 
2/ (X) - [X(l - XT 1/2+£ , 

y(C)-C 1+2£ (i-C) 2£ , 

and 

( r a a, r\ i ^ (l-x)(i-Q , x(i-C) 

P(l-P) a(l-a) 61(1-0) 7(1-7) 

This contribution is finite and yields 

-rtv/^^W . ,< 7 Z- ( 4 ° + 32JV + 8Ar2 + Wa ) [7 .056 1 9233] 
J P (P + ?7 ) »7 3 (4tt) 7 648 

The p — case is given by 

d ' A - /'_ = ;;,„' -,/,;. • — — m +4, . 



U- / S 6 (0) _ x4 w 4 TV (40 + 327V + 87V 2 + TV 3 ) | / M 2 ex P 7 ^ 

648 \ rj 2 

yi(a)y{(3)y{l)yi{e)yi{4>)y{x) f 1 



(p2 + ^2)2 ^2 (4^)6 64g V y y ^ 2 y 



^ W ,7/W ARAM ^IM^ A 



[y(a,p\7,M,x)] 1+4£ 7 P b 2 +^ 2 ) 5 



where 



y(a, P, 7, 0, <j>, x) =4>0 + 



yf(a) = [a(l - a)]" 1 / 2 " 6 , 
= (1 - e)- 1 ^ , 

= 1/2+£ (l - 0) 2e , 

0(1 -fl) , x(W) , (i-x)(i-0) 



a(l-a) 0(1-/3) 7(1-7) 



Integrating one obtains 



S 6 (0) xi u 4 N (40 + 327V + 8iV 2 + N 3 ) 



6 N / 7 2 , 2A2 = -^7T^7- ^ qiO.21524] 

' (p 2 + ?? 2 ) 2 r) 3 (47r) 7 648 
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25 



with 



*8 = 1-£(1-X) + £(1-X)- 



0(i-e) 

' a(l-a) 



Integrating one obtains 



<5 4 /V 



(1-6)[1- </>(!- 6)] 



S 8 (p) r4 u 4 A (44 + 327V + 57V 2 ) r 
- 5 — — — — — ^-[4.31098] 



p (p 2 + rj 2 ) 2 rf (An) 
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The p = case is given by 



(A108) 



(A109) 



<5 4 A 



s 8 (Q) 

p (p 2 +?y 2 ) 2 



= 



TV 44 + 32A + 5N 



K AT 2 \ 



v 3 ( 47r )15/2 
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-r(l+4e)r(l/2 + e) 



M 2 exp7 B 



I dad(3d£d9d<j)dx 
'a 



k(a)k(/3)k(^)k(9)k((j))k(x)k((f>, fl)A, 



-3/2-3e 



(A110) 



with the same notation as used in Eq. (A97). Integrating Eq. (A110) one obtains 



S 4 N 



E 8 (0) 



TV (44 + 32A + 5A 2 ) 



[6.02057] 



/p (p 2 + t 1 2 ) 2 yy 3 (4tt) 7 324 

We have another five-loop contribution given by the twentieth graph of Fig. 3 whose contribution is 



(Alll) 



V / ^ifl = -^^g + f [1.4803] 



p (p 2 +i] 2 ) 2 



7] 3 (18) 2 (4tt) 



(A112) 



where M. = 3 accounts for the three possible ways of inserting one setting sun within another graph of the same type. 
The first term on the LHS of Eq. ( |A112j ) is 



-5 4 /V 



= 8' 



s 9 (p) 

p (p 2 +i] 2 ) 2 



4 A(A + 2) 
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1 [E^-E^O)] 



pkq 



(p 2 + rf) 2 (k 2 + if) (q 2 + ri 2 ) 2 [(p + k + q) 2 + rf 



u 4 A(A + 2) 2 r(3/2 + 5e) /M 2 exp 7£ ; 



rf 



2c 



(18) 2 (4tt) 15 /2 V v 2 
_ 3 x'(a,P,j,9,cf),x) 
2 x(a,f3,"j,6,<j),x) 



5e ,1 



a ara aoaaa x(a)x(P)x(-y)x(9)x((p)x(x) 

o dad ^ ded ^ K^, 7lM ,xF^ 



5 x'(a,(3,-f,9,(j),x) 
2 x(a,0,j,e,(p,x) 



4 u 4 N(N + 2) 2 [n x 10- 
t? 3 (18) 2 (8tt) 2 



i +10 In f — ) - 5.9258 
e V ?7 



(A113) 



Note that £i(<z) and £i(0) have, except for the labeling of momenta, the same form as Eqs. (AE) and (A15) with 
m = j = h = 1. The x functions are given by 



x(a) = [a(l -a)]-^ 2 - e 



(A114) 



^/3) = (l-^)- 1 / 2+e [/3(l-/3)]- 2£ , 



(A115) 



z( 7 ) = (1 - 7) 2 ' 



(A116) 



26 




27 




28 



= [1-0(1- 



(A139) 



A 



g 2 e(i-x) 
[i-0(i-5)P 



flg(i-x) 

- - 0(1 - 61) 



+ 0£(i-x), 



(A140) 



with 



and 



e 



i-x-g(i-x) 

0(1-/3) 



5(1-5) 



0(1 -5)[1- 0(1-0)] 



x + $£(i-x) , 



2(9 



1-0 I-, 



[1-0(1-5)]' 



0(1-5) 



,(1 - a)0(l - 5)[1 - 0(1 - 5)] 0(1 - 5)[1 - 0(1 - 



Integrating, one obtains 



(A141) 



(A142) 



(A143) 



S 4 N 



Sn(p) 

p (p 2 + T] 2 ) 2 



= s 4 



N (44 + 32A + 5N 2 ) 



rf (47r) 7 
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[6.12476] 



The p = case is given by 
Sn(0) 



S 4 N 



p {p 2 +r/ 2 ) 2 



= -S 4 



i 4 N (44 + 32iV + 5A^ 2 



v 3 (4^)15/2 



324 



-r(l +4e)r(l/2 + e) 



./0 



— l+4e 



with the same notation as used in Eq. ( A131 ). Integrating Eq. ( A145 ) one obtains 



(A144) 



M 2 expjE 
rf 



(A145) 



S 4 N 



Sii(O) 
(p 2 + rf 2 ) 2 



N (44 + 32A + 5N 2 ) 



rf (47r) 7 



324 



[8.50217] 
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TABLE I: Comparision of the results for ci as obtained from different methods (see text) and at different orders of approxi- 
mation. 
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1/N NLO 
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FIG. 1: Vertices (top) and diagrams contributing to the self-energy up to order 8 4 . 





FIG. 3: All diagrams contributing to the two-point function ((j) 2 )s, up to order S 4 , at the critical point. Again, the black dot 
represents here only the Sif insertions. 



